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ABSTRACT 

Diffusing species into some host metals in semi-infinite cylindrical solids and key thermodynamic properties were examined. The 
governing Bessel's differential equation of order zero and time dependent first order differential equation were obtained. Analytical 
solutions showed that stable concentration of the system and partition function was achieved as the temperature increases. 
However, the thermodynamic function of Helmholtz free energy is enhanced as the radius of the diffusing species into the host 
metals at a fractional exponent of 0.6 is increased and the reverse was observed with the entropy of the system while the internal 
energy and specific heat are independent of the diffusing species parameters. A mathematical model was also proposed to calculate 
the diffusion coefficients of the diffusing species into the host metals by adjustment of the constant in the model. 
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1. INTRODUCTION 


Diffusion is regarded as the stepwise movement of atoms from lattice site to lattice site. It can also be viewed as the process of 


mixing which involves the migration of atoms from region of higher concentration to that of lower concentration or the shifting of 
atoms and molecules to new areas within a given material resulting in the uniformity of composition as a result of thermal agitation 
(Rajput2013). Diffusion is fundamentally statistical physics in nature and the process irreversible hence increases entropy. It is also 
synonymous with phase change and movement is short-range. Diffusion in solids is a slow process and its rate increases with 
increase in temperature which is vital in metallurgical processes. As a result of the applications of diffusion which include phase 
changes, metal bonding, oxidation of metals, doping of semiconductors, recovery and recrystalization and a host of researches are 
abounded. Hung et al. (2000), in the study of self-diffusion in metals using the fourth order moment approximation, deduced the 
activation energy and pre-exponential values for some metals at high and low temperature regimes which are in good agreement 
with experimental results. Martinez-Vera et al. (2015), examined the diffusion co-efficient estimation in shrinking solids and 
presented an experimental results of the drying kinetics for a solid that shrinks in volume during the drying process depending on 
the drying temperature. Povstenko et al. (2017) tackled fractional diffusion in a solid with mass absorption wherein, a model was 
used to illustrate far reaching results. Carballido et al. (2011), proposed a turbulent diffusion of large solids in a protoplanetary disc 
with two models to predict the same radial turbulent diffusion coefficient for small particles. Duda et al. (2007), in the study of 
modeling of coupled deformation- diffusion — damages in elastic solids, used numerical algorithm to solve the discrete system of 
equations for the application of the hydrogen embrittlement. In the study of volumes of activation for diffusion in solids, Reyes 
(1958), critically assessed the activation volumes of some simple models of lattice defects and compared with experimental 
measurements. The comparison was satisfactory in many cases. A general mathematical approach to the rigorous treatment of 
experimental data on non stationary state diffusion in solids of complex shape was developed by Beeker (1959) in the study of 
diffusion in solids of arbitrary shape with application to the drying of the wheat kernel. A new theoretical approach was established 
by Nguyen (2018) to define transport coefficients of charge and mass transport in porous material directly from impedance data ina 
study of new approach to determine gas diffusion coefficient in porous solids by Electrochemical Impedance Spectroscopy (EIS): 
Application for NH3 and CO2 absorption Zircomia and Zeolite type 5A. In the study of measurement of solids diffusion coefficient by 
a temperature programmed method, Kapoor and Oyama, (2011), presented a method for determining diffusivities in solids where 
the diffusing species absorbs or reacts at the external surfaces and where the diffusivity does not vary appreciably with 
concentration. The method requires measuring the flux of the diffusive species out of the solids under the influence of temperature 
program. A model is proposed based on isothermal Fickian diffusion. Korblein et al (1985), examined the diffusion process in solid Li 
-Mg and Li - Ag alloys and the spin-lattice relaxation of ®Li and opined that Ag is found to diffuse faster than Li. Experimental results 
displaying the quantum diffusion in different systems was the plank of the study of quantum diffusion in solids by Kagan (1992). Ina 
study of sucrose diffusion in aqueous solution, Price et al (2016), reported measured sucrose diffusion coefficient using Stokes- 
Einstein equation and claimed it fits well with experimental data. A recent study of transfer phenomena in fluid and heat flows VIII by 
Kolisnychenko (2019) highlighted a compilation of theoretical and practical aspects of modeling and numerical investigation of 
diffusion convection and MHD mixed convection flows. The main aim of this study is to consider temperature and concentration 
dependent systems of diffusion coefficient and highlight some of its thermodynamic properties in semi-infinite cylindrical solids. 


2. FORMALISM 


The governing partial differential equation of a long circular cylinder in which diffusion is everywhere radial, concentration C is 
therefore a function of radius r and time t and the model becomes 


cc 10 [ <) 
= rD (1) 
Oot ror or 


with the boundary conditions 


CA anec, 
C(2,t) =C, (2) 
C(r,0) =C, 
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The diffusion coefficient D (concentration- dependent system) is given by 


= AC 
D=D,e% (3) 
where D, is the diffusion coefficient for the surface concentration of the solids and B, is a constant. 


The diffusion coefficient in a temperature-dependent systems obey the Arrhenius-type equation 
E, 


D= Dye * (4) 


where E.., R, T and D, are respectively the activation energy, universal gas constant, temperature and pre-exponential factor. 


A generalized expression for the activation energy according to Atkins and Paula (2006) is given as 


dInD 
(5) 


bak 
dT 
which incorporate the quantum mechanical tunneling role in the diffusion of the reacting metals. 


According to Price et al (2016), the Fractional Stokes-Einstein relation for diffusion constant rate is given as 


kT 
6zur* 


D= (6) 


where € <1 and ranges from 0.6 —0.9, su is viscosity and k, is Boltzmann constant. 


To determine the activation energy of the system, equation (6) is put into equation (5) and the result is given as 


E, =RT=N,k,T (7) 


where NV, is Avogadro's constant. 
The quantity of diffusing substance Q, which diffuses through unit length of a cylinder in time t is given by Crank (1975) as 


= 2nDdC, =) 


In(2) 


If Q, is measured in a concentration - dependent system, the mean value of the diffusion coefficient from equation (8) is given as 


Cy 
| Ddc 
q 
(C, = Ci) . 


Method of solution 


To solve equation (1), separating the variables techniques is assumed, which is of the form 
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C(r,t)= R(T) (10) 


Equation (10) is put into equation (1) and the resulting expression is equated to a constant — Bg , two equations result. They are 


T'(t)+ B’DT(t) =0 (11) 
1 
R'(r)+—R(r)+ RIB’ =0 (12) 
r 
Using the method of variable separable for equation (11) and the Frobenius method of the form R(r) = Lar for equation 


n=0 
(12) and substituting the resulting equations into equation (10), the result is given as 


2.2: 4_4 23 405 
C(r,t)= air sae Ms a(n? et +] (-#) (13) 


4 64 9 225 


where the constant of equation (11) has been absorbed into Gy and d,. 


2 
CG. "O66; 


Imposing the boundary conditions and solving the resulting simultaneous equations, the values of GQ) = and 


a, 
7 C,a, — aC, 


Bit, — Aa; 


where @, -|[0&s Fe el (e-#") 

a =| Oo Bo +]fer) 
fet 8 Jer" 
“| 


[282 BE ‘| (.-#) 


To determine the concentration from the concentration-dependent diffusion coefficient, equation (3) is equated to equation (6) and 
the expression takes the form 


1 k,T 
C=—ln —— (14) 

B, D, 67ur 
Following Atkins and Paula (2006), equation (6) is put into equation (3), the expression is put into equation (4) and finally into 
equation (5), the diffusion coefficient takes the form 


2 k T 2 
D= Exp - In B . (15) 
2 \D,6aur’) 4 


where the constant of integration is assumed unity 
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3. THERMODYNAMIC FUNCTIONS 


It has been stated earlier in the introduction that diffusion is statistical physics in nature, therefore to determine the thermodynamic 


properties of the diffusing species into host metals cannot be overemphasized. Any physical system with thermodynamic properties, 
needed the introduction of the concept of partition function given by the expression 


“ E 
2,0) =¥Ev-( 7 ) (16) 
B 


n=0 


where E’, is energy. 


Equating equations (4) and (6), the expression takes the form 


T 
E, =-RT In| —“2? _ (17) 
D,6anr* 


Assuming the activation energy to be the energy of the partition function, equation (16) takes the form 


D,6anr* 


n=l 


é T 
Z,(T) = > ExpN, of 5) (18) 


Thermodynamic functions such as Helmholtz Free Energy F(T), Entropy S(T). Internal Energy U(T) and Specific Heat C, (7) are 


obtained from the Partition function as follows 


F(T)= -EInZ(D) (19) 
OF (T) 

S(T) =-k, ——— 20 
(T) Bar (20) 
U(T)= =, Z(T) (21) 

oT 
C,(T) =k, ao (22) 


Values of constants and other parameters used 


k, =1.38065x10 J /k 
N, =6.02214x10” mol” 
R=8.31447J /kmol 


C, =lkg/mol 
C, = 2kg/mol 
é€=0.6 
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Table 1: Diffusion Data (Rajput, 2003; Hilderbrand and Lamoreaux, 1976) 


—_ Viscosities Pre-exponential Factor | Activation nee - 
Diffusing Species n(Ns | m?) Host Metals D, (m? /s) wets je i 
Copper (Cu) 0.44 Cu 7.8x10> 211 4.2x10°"° 
Zinc (Zn) 0.18 Cu 2.4x10~ 189 4.0x107'8 
Aluminium(Al) 0.09 Al 9.3x10~ 144 4.2x10-4 
Copper (Cu) 0.44 Al 1.2x10~ 131 4.0x107'* 
Copper (Cu) 0.44 Nickel(Ni) 2.7x107 256 1.3x10~ 


Table 2: Calculated values of diffusion parameters 


Diffusing Species-Host Concentration Constant Diffusion coefficient 
Metals C(kg/s) PB, (mol/kg) D(m* /s) 

Cu-Cu -38.6433 0.850 4.2354. x107” 
Zn-Cu -36.5709 0.803 4.23186 x10 8 
Al-Al -38.1377 0.588 4.19458 x10~'4 
Cu-Ni -37.5828 1.061 1.29974 x10 ~~ 
Cu-Al -38.4610 0.551 4.06762 x10~'4 


1 380 
1 385 
1 390 


1 3995 


Figure 1: Dependence of Concentration on radius of diffusing entity with Temperature varying for Cu—Cu diffusion 
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Figure 4: Dependence of Concentration on radius of diffusing entity with Temperature varying for Cu — Al diffusion 
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Figure 5: Dependence of Concentration on radius of diffusing entity with Temperature varying for Cu — Ni diffusion 
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Figure 6: Dependence of Helmholtz Free Energy on Temperature with radius of diffusing entity varying for Cu-Ni diffusion 
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Figure 7: Dependence of Entropy on Temperature with radius of diffusing entity varying for Cu—Ni diffusion 
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Figure 8: Dependence of Concentration on radius of diffusing entity with Temperature varying for Cu—Cu diffusion 


4. DISCUSSION 


Equation (14) and Figures (1-5), depict the behaviour of concentration of diffusing species on temperature. It is shown that increase 
in temperature of the diffusing species result in a corresponding increase in the concentration. This observation indicates that 
increase in temperature, provide higher diffusion coefficients because the atoms have higher thermal energies and therefore higher 
certainty of being activated over the energy barrier between the atoms. This observation is in agreement with the work of Rajput 
(2003). Tables 1 and 2 also show that the values of diffusion coefficient obtained from our study fits well with that of experimental 
data. Equation (18) gives an expression for the partition function which described the thermodynamic properties of the diffusing 
species and host metals of the system. The equation showed that increase the radius of diffusing species at a fractional exponent of 
0.6 within the temperature range of 773k — 1573k, the partition function appears to be stable in Cu-Cu diffusion, Zn-Cu diffusion, Al- 
Al diffusion, Cu-Al diffusion and Cu-Ni diffusion. The stability is associated with the high Avogadro’s number. Figure (6) was 
modeled from equation (19). It is observed that increase in the radius of the diffusing specie at 0.6 exponent enhanced the 
Helmholtz free energy of the Cu-Cu diffusion. The same condition is observed with the Zn-Cu diffusion, Al-Al diffusion, Cu-Al 
diffusion and Cu-Ni diffusion. Figure (7) is the graph of equation (20). It depicts the dependence of the entropy on the temperature 
of the diffusing specie. An increase in the radius of the diffusing specie diminished the rate of disorderliness of the system. This 
observation is the same with the other diffusing species and host metals. Equations (21) and (22) respectively give an expression of 
the internal energy and specific heat at constant volume of the systems under consideration. The equations show that these 
thermodynamic functions are dependent on the temperature but independent of the pre-exponential factors, viscosities of the 
diffusing species and radius of the diffusing species. The simple explanation is that the internal energy and specific heat are not 
affected by the mentioned properties of the diffusing species. Equation (13) and graph (8) explained the relationship between the 
concentration of the systems and radius of the diffusing specie. A stable concentration of the system is achieved as the temperature 
increases. Although, the constant factor @ characterize the different diffusing species and host metals in magnitude, the overall 


observation as the temperature increases for all the solid metals under consideration is stability after the initial increase in the 
concentration of the systems. A mathematical model for determining the constant diffusion coefficients for the solid metals under 


consideration is proposed. Analysis showed that the constant B, varies from the different diffusion species and host metals shown 
on Table 2. The value of the concentration for a given diffusion specie is obtained from equation (14), substitute into equation (3) 
with the appropriate constant Bus the diffusion coefficient for the diffusing specie and host metal under consideration can be 


obtained. 


5. CONCLUSION 

The study examined an analytical method in the consideration of the diffusing species into host metals in cylindrical coordinate. The 
plank of the study is to lay credence to the assertion that the fractional -Stokes-Einstein equation is a better approximation to the 
Stokes-Einstein equation in the determination of diffusion coefficients of solids. Results obtained in the study and in particular the 
model proposed is in reasonable agreement with most of the literatures cited and fits experimental data. 
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